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EFFECT  OF  UNIAXIAL  STRESS  ON  THE  RAMAN  SPECTRA  OF 

GRAPHITE  FIBERS 

H.  Sakata^,  G.  Dresselhaus  and  M.S.  Dresselhaus 
Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139 

M.  Endo 

Shinshu  University,  Nagano-shi  380,  Japan 


Abstract 

Measurements  of  the  effect  of  uniaxial  stress  on  the  frequency  of  the  it  as  0  opti¬ 
cal  phonons  in  heat-treated  benzene-derived  graphite  fibers  using  first-order  Raman 
scattering  have  been  performed.  Application  of  uniaxial  stress  along  the  fiber  axis  is 
found  to  cause  polarization-dependent  splittings  and  shifts  of  the  Raman  peaks  of  the 
k  as  0  doubly  degenerate  optical  phonons.  From  these  observed  splittings  and  shifts, 
experimental  values  for  the  phenomenological  coefficients  which  describe  the  changes 
in  the  elastic  constant  of  these  phonons  with  strain  are  determined.  It  is  concluded 
that  Raman  spectroscopy  can  be  used  to  characterize  the  local  stress  or  strain 
conditions  of  graphite  fibers  nondestructive^. 

Introduction 

Laser  Raman  spectroscopy  is  a  powerful  nondestructive  technique  for  monitoring  the 
local  stress  variations  near  the  surface  of  materials.  The  Raman  microprobe  used  in 
this  experiment  provides  ~  2jxm  spatial  resolution  within  the  optical  skin  depth.  This 
spatial  resolution  is  to  be  contrasted  with  more  conventional  techniques  such  as  x-ray 
analysis  which  typically  have  a  spatial  resolution  of  ~  1mm.  Stress  measurements 
in  semiconductors  such  as  Si[l],  Ge[2]  and  GaAs[2]  by  Raman  spectroscopy  have  been 
previously  reported.  Of  particular  interest  for  Raman  spectroscopic  studies  of  graphite 
is  the  anisotropy  in  the  stress  and  strain  parallel  and  perpendicular  to  the  basal  planes, 
so  that  polarization  studies  provide  unique  information  on  the  elastic  properties. 

Graphite  fibers  are  currently  of  interest  to  materials  applications  requiring  high 
strength  to  weight  ratio.  However,  no  stress  measurements  in  graphite  fibers  by  Raman 
spectroscopy  have  been  reported,  though  Raman  measurements  of  the  in-plane  optical 
lattice  mode  have  been  shown  to  provide  a  sensitive  method  for  the  characterization 
of  the  structural  perfection  of  carbon[3]  and  carbon  fibers[4]  [5].  Since  the  benzene- 
derived  graphite  fibers[6]  are  typically  between  10  and  20/im  in  diameter,  Raman 
spectroscopy  can  be  used  to  examine  the  stress  variation  within  the  optical  skin  depth 
(~  600A  for  light  scattering  at  4880A). 

In  this  work,  heat-treated  benzene-derived  graphite  fibers  are  characterized  by 
Raman  spectroscopy  as  a  function  of  applied  uniaxial  tensile  stress.  In  the  absence 

^Permanent  address:  Mechanical  Engineering  Research  Laboratory,  Hitachi  Ltd.,  Tsuchiura-shi, 
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of  uniaxial  stress,  the  first-order  Stokes  Raman  spectrum  of  graphite  fibers  exhibits  a 
single  peak  which  corresponds  to  the  few  0  vibrations  of  the  doubly  degenerate  opti¬ 
cal  phonons  (Eigz).  The  application  of  uniaxial  stress  causes  polarization-dependent 
splittings  and  shifts  which  are  observed  to  be  linear  in  the  applied  stress  in  consequence 
of  the  change  of  the  elastic  constants  with  strain.  From  these  observed  splittings  and 
shifts,  experimental  values  for  the  phenomenological  coefficients  .4  and  B  which  de¬ 
scribe  the  changes  in  the  elastic  constant  of  the  k  %  0  optical  phonons  with  strain  are 
determined. 

Experimental  Procedure 

The  carbon  fibers  used  in  the  present  investigation  were  prepared  on  a  substrate  by 
pyrolyzing  a  mixture  of  benzene  and  hydrogen  at  a  temperature  of  ~  1100°C.  After 
the  fiber  growth,  graphitization  of  the  fibers  is  accomplished  by  heat  treatment  at 
a  temperature  of  2900°C  using  a  carbon-resistance  furnace  in  a  high-purity  argon 
atmosphere.  The  resulting  fibers  consist  of  faceted  concentric  graphite  layers  around 
the  fiber  axis,  such  as  the  annular  rings  of  a  tree  [6].  The  fiber  diameter  used  in  this 
work  was  7.6pm  and  the  fiber  length  was  15  mm.  As  shown  in  Fig.  1,  single  fibers 
were  bonded  to  the  sample  stage  of  the  Raman  microprobe  at  one  end  and  to  a  nylon 
string  which  was  connected  to  a  pan  containing  weights  at  the  other  end.  Uniaxial 
tensile  stress  was  applied  to  the  fibers  by  putting  weights  on  the  pan  and  thus  the 
stress  direction  was  coincident  with  the  fiber  axis. 

The  Raman  scattering  measurements  were  made  at  room  temperature  in  the  backscat- 
tering  configuration  using  a  Raman  microprobe  [7]  with  488  nm  Ar+  laser  excitation, 
a  40A"  objective  lens  (about  1pm  spot  size)  and  a  diode  array  detector.  The  Raman 
spectra  were  taken  for  polarizations  such  that  the  electric  vectors  of  the  incident  beam 
were  either  parallel  or  perpendicular  to  the  stress  direction  of  the  fibers.  Spectra  were 
also  taken  for  more  than  8  points  at  random  along  the  fiber  axis  at  each  stress  level. 
The  values  of  the  peak  position  and  half-width  in  the  Raman  spectra  were  determined 
by  a  Lorentzian  fit  to  the  experimental  points. 

Experimental  Results 

The  doubly  degenerate  Raman  mode  at  about  1580cm-1  which  corresponds  to  the 
Raman-allowed  Ei92  mode  in  HOPG  was  used  for  the  Raman  characterization.  The 
disorder-induced  line  at  about  1360cm-1  could  barely  be  detected  in  the  Raman  spec¬ 
tra  of  the  fibers  used,  indicating  a  nearly  complete  graphitization  of  the  fibers  and  a 
nearly  full  establishment  of  three-dimensional  graphite  ordering[5]. 

Figure  2  shows  typical  Raman  spectra  of  the  fibers  for  the  uniaxial  tensile  stress 
<7  =  0  and  430  MPa.  The  stress  a  is  applied  along  the  fiber  axis,  while  the  electric 
vector  E  of  the  incident  beam  is  polarized  either  parallel  or  perpendicular  to  the 
stress  direction,  that  is.  a\\E  and  <f±E,  respectively.  For  the  fiber  batch  used  in 
these  experiments,  the  Raman  shift  of  the  peak  in  the  absence  of  stress  occurs  at 
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wo  =  1583.2cm-1  and  its  half-width  To  at  half  maximum  intensity  is  10.3  cm-1.  From 
Fig.  2,  it  can  be  seen  that  the  application  of  uniaxial  stress  causes  the  Raman  shift  for 
the  peak  to  split  into  two  components  which  are  observed  for  the  two  polarizations. 

Figure  3  shows  a  plot  of  the  Raman  shift  for  the  peak  near  1580  cm-1  as  a  function 
of  the  uniaxial  tensile  stress  a.  Since  spectra  were  taken  for  more  than  8  points  at 
random  along  the  fiber  axis  at  each  stress  level,  mean  values  and  variation  bars  are 
shown  in  Fig.  3.  The  Raman  shift  of  HOPG  which  was  obtained  in  this  experiment 
is  also  presented  in  this  figure  as  a  reference  and  it  can  be  seen  that  the  Raman  shift 
of  the  fibers  in  the  absence  of  stress  is  within  1  cm-1  of  that  of  HOPG  (w  =  1583.6 
cm-1).  From  our  data  in  Fig.  3,  it  follows  that  the  Raman  shift  decreases  linearly  with 
increasing  a  in  both  cases  {3\\E  and  a±E)  and  the  decrease  in  the  case  of  5\\E  is  about 
3.5  times  faster  than  that  of  oLE.  The  different  behavior  for  the  two  polarizations 
corresponds  to  the  splitting  of  the  double  degeneracy  of  the  E^g?  mode  under  uniaxial 
stress. 

The  variation  in  the  half-width  T  at  half  maximum  intensity  of  the  1580cm-1 
lines  as  a  function  of  a  for  the  two  polarization  directions  is  shown  in  Fig.  4.  The 
half-width  To  of  the  fibers  in  the  absence  of  stress  is  about  10%  larger  than  that  of 
HOPG  (T  =  9.1  cm-1).  From  Fig.  4,  it  can  be  seen  that  the  half-width  T  increases 
linearly  with  increasing  a  for  both  polarizations  and  the  increase  in  the  case  of  5\\E 
is  about  3  times  faster  than  that  of  dL E.  These  half-width  variations  correspond  to 
the  inhomogeneous  variation  of  the  local  stress  or  strain  distribution. 

Figure  5  shows  a  plot  of  the  shift  Aw  =  wo  -  w<,,  where  w0  and  w*  are  the  Raman 
shifts  in  the  absence  and  presence  of  stress,  respectively,  as  a  function  of  a.  The  data 
for  diamond  (ct||[100],(7_L.E)[8]  and  semiconducting  materials  such  as  Si [  1  ] ,  Ge[2]  and 
GaAs[3]  (<7||[lll],<?||i?)  are  also  presented  in  Fig.  5.  From  Fig.  5.  it  follows  that  Aw  for 
the  fibers  increases  linearly  with  increasing  tr  for  both  polarizations  and  the  increase 
of  Aw  in  the  case  of  S\\E  is  about  3.5  times  larger  than  that  of  alE.  Thus  the  stress 
<t  is  related  to  Aw  by  the  following  equations: 

tj(MPa)  =  168  •  Aw(cm-1 );  cr\\E  ( 1  < 

a(MPa)  =  588  ■  Aw(cm-1);  <r±E  (2 

These  experimental  measurements  should  be  useful  in  engineering  applications  since 
the  local  stress  can  be  estimated  nondestructively  by  measurement  of  Aw.  It  can  also 
be  seen  that  the  increase  of  Aw  for  the  benzene-derived  graphite  fibers  is  comparable 
to  that  for  other  materials  shown  in  Fig.  5.  We  have  been  unable  to  obtain  this 
measurement  on  crystalline  graphite. 

The  relationship  between  Aw  and  the  strain  e  which  is  calculated  from  the  elastic 
modulus  and  the  applied  stress  is  shown  in  Fig.  6,  where  it  is  found  that  the  strain 
dependence  of  Aw  for  the  graphite  fibers  is  somewhat  larger  but  comparable  to  that 
for  typical  semiconducting  materials. 

From  this  Raman  scattering  study  of  the  characterization  of  heat-treated  benzene 
derived  graphite  fibers  under  stress,  it  is  concluded  that  Raman  spectroscopy  can  be 


used  to  provide  complementary  information  on  the  local  stress  conditions  of  these  j 

fibers.  | 


Theoretical  Considerations 

1  Stress— Strain  Relations  for  Graphite  Fibers 

For  hexagonal  crystals,  symmetry  reduces  the  generalization  of  Hooke’s  law  which 
relates  an  arbitrary  stress  to  an  arbitrary  strain  as  follows 
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The  carbon  fibers  used  in  this  experiment  exhibit  a  unique  structural  morphol¬ 
ogy,  since  the  carbon  was  deposited  in  a  cylindrically  symmetric  fashion  leading  to  a 
smooth  uniform  lamellar  deposit  in  which  the  graphitic  basal  planes  had  a  high  degree 
of  preferred  orientation,  lying  nearly  parallel  to  the  fiber  axis[6].  Hence,  taking  the 
coordinate  system  shown  in  Figure  7,  for  a  uniaxial  stress.  <r,  applied  along  the  fiber 
axis, 

axx  =  a  (4) 

Gyy  =  °zz  ~  Gy z  =  G  zx  —  GXy  =  0  (5) 

Also,  for  a  uniaxial  stress  a  there  are  no  shear  components  of  strain,  so 

iyz  =  izx  =  ixy  =  0  (6) 

Therefore,  the  stress-strain  relations  for  graphite  fibers  are 

i xx  =  S\\G 

iyy  =  S\2G 

fzz  =  SnCT  (7) 


2  Strain-Dependent  Raman  Frequency  in  Graphite 


In  the  presence  of  strain  the  dynamical  equations  which  describe  phonon  modes  in  the 
solid  have  the  form[l],[2],(8] 

fniii  = (8) 
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where  u;  is  the  ittl  component  of  the  relative  displacement  of  the  two  atoms  in  the  unit 
cell,  m  is  the  reduced  mass  of  the  two  atoms,  and  K°  =  mwj  is  the  effective  spring 
constant  of  the  phonon  modes  in  the  absence  of  strain.  The  double-dotted  quantities 
indicate  second  derivatives.  The  quantity 

dlv' 

~fc~€kl  =  K'ijklekl  =  K'ijlk€tk  (10) 

is  the  change  in  the  spring  constant  due  to  the  applied  strain  i,j,k  and  l  denote 
crystallographic  axes. 

From  Eq.  9  the  phonon  frequencies  for  the  strained  lattice  can  be  deduced.  For 
pristine  graphite  with  4  atoms/unit  cell  a  12  X  12  secular  determinant  is  obtained  for 
the  3  acoustic  and  9  optical  modes.  The  hexagonal  symmetry  of  the  lattice  results  in  a 
factorization  of  the  secular  determinant  at  k  =  0  and  the  result  for  the  Raman  active 
E232  mode  is  given  by 

A(€xx  F  £yy)  ~  A  B((.Xx  ~  tyy  F  2ifry) 

B{(xx  ~~  tyy  ~  2 i^xy)  4(frr  F  Cyy)  ~  ^ 


where 


A  =  w*  -  u>q. 


•2  (12) 

The  parameter  is  the  strain-dependent  Raman  frequency  and  wo  is  the  frequency 


in  the  absence  of  strain.  The  above  can  be  approximated  as 
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Since  the  shear  components  of  strain  vanish,  the  secular  equation  reduces  to 

=  0. 


A(eXx  F  Cyy  )  A  B(^€xx  fyy ) 
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(14) 


Hence  there  are  two  solutions: 

Al  =  (4  F  B)(xx  F  (4  —  B)tyy 
A2  =  (4  ~  B)(XX  F(AF  B)(yy. 


(15) 


Substituting  the  stress-strain  relations  (Eq.  7)  into  Eq.  15,  the  stress-dependent 
Raman  frequencies  become 


w„i  =  wo  F  - — {(4  +  B)S ii  F  (4  -  .BJSn} 
2wo 


w<,2  =  ^0  F  - — {(4  -  B)S  11  F  (.4  F  B)5j2}. 
Zwq 


(16) 


The  shifts  in  the  Raman  frequency  due  to  strain,  Awi  and  AW2,  are 


Au>i  =  wo  —  w<,i  —  —  - — {(/l  F  5)5jj  F  ( 4  -  B)S\2} 

Zwo 
a 


Aw2  =  wo  —  w„2  =  -- — {(4  —  B)S\\  F  (4  F  B)S\2}. 

ZWq 


(1~) 
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Using  the  values  for  the  elastic  compliance  constants  in  graphite,  which  have  been 
determined(9]  and  the  experimental  data  for  the  relationship  between  the  applied  uni¬ 
axial  stress  a  and  the  shift  Asji  or  Au2,  the  phenomenological  coefficients  A  and  B 
which  describe  the  changes  in  the  elastic  constant  of  the  k  =  0  optical  phonons  with 
strain  are  determined: 

A  =  -1.44  X  10‘cm-2 

B  =  ±5.80  x  106cm-2  (18) 

From  the  above,  it  is  clear  that  the  Raman  effect  can  be  used  to  characterize  the 

strain  induced  in  a  carbon  fiber  which  is  subjected  to  an  unknown  stress  load,  a, 

either  externally  or  internally.  If  polarized  light  (E  ||  a  and  E  ±  8)  is  used  for  the 
Raman  characterization,  the  shift  of  the  frequency  of  the  E2S2  Raman  active  mode  will 
give  information  on  the  magnitude  of  the  stress  (or  the  strain)  and  its  homogeneity 
along  the  fiber  axis  on  a  ±2/im  resolution  scale,  using  Eqs.  17  and  18.  From  the 
linewidth,  additional  information  is  provided  on  the  uniformity  of  the  local  strains. 
If  unpolarized  light  is  used,  the  shift  of  the  E2ff2  Raman-active  mode  A u  will  be  a 
linear  combination  of  Awi  and  Au>2,  which  for  completely  random  polarization  becomes 
Au  =  (Au;i+Au;2)/2.  In  this  case,  the  linewidth  will  mostly  reflect  the  superposition  of 
the  contributions  of  the  two  polarization  contributions  which  may  not  be  well  resolved. 

Conclusions 

1.  Raman  spectroscopy  can  be  used  to  provide  complementary  information  on  the 
local  stress  conditions  of  graphite  fibers.  The  local  stress  can  be  estimated  non- 
destructively  by  measurement  of  Au. 

2.  The  stress  dependence  of  the  first-order  Raman  spectra  associated  with  the  k  x  0 
optical  phonons  in  graphite  fibers  has  enabled  us  to  determine  values  for  the 
phenomenological  coefficients  which  describe  the  changes  in  the  spring  constant 
of  these  optical  phonons  with  strain. 
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Figure  Captions 

Figure  1:  Schematic  drawing  of  the  experimental  setup  for  Raman  scattering. 

Figure  2:  Line  shape  of  the  optical  phonon  E2S2  of  the  fibers  for  the  uniaxial  stress 
<7  =  0  and  430  MPa.  The  solid  line  is  the  best  Lorentzian  fit  to  the  experimental 
points.  The  dashed  line  indicates  the  position  of  the  Raman  shift  a>o  in  the  absence 
of  stress.  The  arrows  indicate  the  positions  of  the  Raman  shift  for  the  peak  in  the 
presence  of  stress.  In  the  case  of  a  =  430MPa,  the  Raman  shift  of  the  peak  occurs  at 
u ;  =  1580.6  and  1582.4  cm"1  and  its  HWHM  linewidth  P  is  11.4  and  10.6  cm-1  for 
a\\E  and  a±E,  respectively. 

Figure  3:  Plot  of  the  Raman  shift  for  the  peak  near  1580cm-1  vs.  the  applied  uniaxial 
tensile  stress  a.  The  full  and  open  circles  indicate  the  mean  values  for  the  Raman  shift 
of  the  peak  for  the  graphite  fibers  for  the  two  senses  of  polarization.  The  error  bars 
are  also  given.  For  comparison  the  value  for  HOPG  is  included. 

Figure  4:  Plot  of  the  half-width  T  for  the  peak  near  1580cm-1  vs.  the  applied  uniaxial 
tensile  stress  a.  The  full  and  open  circles  indicate  the  mean  values  for  the  HWHM 
linewidth  for  the  two  senses  of  polarization.  The  error  bars  are  also  given.  Also 
included  is  the  value  for  HOPG. 

Figure  5:  Plot  of  Au  =  uiq  —  vs.  the  applied  uniaxial  tensile  stress  a  for  graphite 
fibers  (G.F.)  including  the  two  senses  of  polarization  of  the  light  and  the  corresponding 
results  for  typical  semiconductors.  The  full  and  open  circles  indicate  the  mean  values 
for  the  Aw  measurements  and  the  error  bars  are  also  given. 

Figure  6:  Plot  of  the  experimental  shift  A u  vs.  the  calculated  strain  e,  which  is  esti¬ 
mated  from  the  elastic  modulus  for  graphite  (including  both  a\\E  and  5LE  polariza¬ 
tions).  For  comparison,  the  corresponding  results  for  Si,  Ge  and  diamond  are  included. 

Figure  7:  Coordinate  system  for  a  crystallite  in  a  graphite  fiber.  The  x  axis  is  taken 
parallel  to  the  fiber  axis  and  is  coincident  with  the  direction  of  the  applied  uniaxial 
tensile  stress  a. 
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